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Abstract 

It is well known that a current driven through a two-dimensional electron gas with Rashba spin- 
orbit coupling induces a spin polarization in the perpendicular direction (Edelstein effect). This 
phenomenon has been extensively studied in the linear response regime, i.e., when the average 
drift velocity of the electrons is a small fraction of the Fermi velocity. Here we investigate the 
phenomenon in the nonlinear regime, meaning that the average drift velocity is comparable to, or 
exceeds the Fermi velocity. This regime is realized when the electric field is very large, or when 
electron-impurity scattering is very weak. The quantum kinetic equation for the density matrix of 
noninteracting electrons is exactly and analytically solvable, reducing to a problem of spin dynamics 
for “unpaired” electrons near the Fermi surface. The crucial parameter is 7 = eELg/Ep, where 
E is the electric field, e is the absolute value of the electron charge, Ep is the Fermi energy, and 
Lg = h/{ 2 ma) is the spin-precession length in the Rashba spin-orbit field with coupling strength 
a. If 7 <C 1 the evolution of the spin is adiabatic, resulting in a spin polarization that grows 
monotonically in time and eventually saturates at the maximum value n{a/vp), where n is the 
electron density and vp is the Fermi velocity. If 7 ;§> 1 the evolution of the spin becomes strongly 
non-adiabatic and the spin polarization is progressively reduced, and eventually suppressed for 
7 —)• 00. We also predict an inverse nonlinear Fdelstein effect, in which an electric current is driven 
by a magnetic field that grows linearly in time. The “conductivities” for the direct and the inverse 
effect satisfy generalized Onsager reciprocity relations, which reduce to the standard ones in the 
linear response regime. 
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I. INTRODUCTION 


The generation of spin polarization by an electric current and, conversely, of an electric 
current by a non-equilibrium spin polarization^^^^ are topics of great interest in spintron- 
ics^^“^®. Both effects have a common origin in the spin-orbit interaction and, in the linear 
response regime, are connected by an Onsager reciprocity relation. Experimentally, current- 
induced spin polarization has been observed in numerous experiments on doped semicon- 
ductors^’^’^^’^^. The inverse effect, known as spin galvanic effect, has been demonstrated in 
semiconductors^^ and very recently in metallic structures®. 

On the theoretical side the spin-polarization effect was theoretically predicted in Refs.^’^ 
in the context of the two-dimensional electron gas with Rashba spin-orbit coupling (Rashba 
2DEG). For this reason the effect is widely known in the literature as Rashba-Edelstein 
effect. Its inverse was studied theoretically in Ref.^ and, more recently, in Ref.^^ (for a more 
complete discussion see Ref.^^. The Edelstein effect also bears a close relationship to the 
theoretically and practically important spin Hall effect.^® Indeed, the spin Hall current in 
the clean Rashba 2DEG arises as a transient in the process of building the Edelstein spin 
polarization. 

All of the theoretical studies mentioned in the previous paragraph were limited to the 
linear response regime - weak electric field or weak spin injection - meaning that the drift ve¬ 
locity of the electrons remains much smaller than the Fermi velocity and the non-equilibrium 
spin polarization is small. There are good reasons for this choice, since this is in practice the 
regime in which virtually all of the experiments have been done. The presence of impurity 
scattering limits the electron drift velocity to values much smaller than the Fermi velocity. 

By contrast, in this paper we present a theoretical study of the Edelstein effect and its 
inverse in a perfectly clean Rashba 2DEG. The absence of impurities allows the electrons to 
be accelerated to high velocities (comparable to the Fermi velocity) and thus to access the 
nonlinear regime. There are several reasons for undertaking this study. First of all, the model 
admits an elegant completely analytical solution, which is no common occurrence in this area 
of research. Second, the solution is very instructive, bringing forth an unexpected connection 
with the classic Landau-Zener-Majorana modeh®“^^ for the anti-crossing of two energy levels. 
In brief, we find that when the drift velocity of the electrons reaches a sufficiently high value 
over a long time (i.e., for weak electric held) the Edelstein spin polarization (normally 
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proportional to the electric field) saturates to a limiting value corresponding to 100 % spin 
polarization of the electrons in an annulus of momentum space comprised between the Fermi 
momenta of the two Rashba bands. However, if the acceleration is very high, the electron 
spins are unable to respond to a rapidly changing spin-orbit field, and the final polarization 
is much smaller than the saturation limit. In addition, we find that, no matter how small 
the electric field is, a Landau-Zener anti crossing always occurs, for sufficiently large times, 
on part of the Fermi surface. States on this part of the Fermi surface can either stay on the 
adiabatic track, or undergo a diabatic crossing, in which case their contribution to the final 
spin polarization is greatly suppressed. 

In spite of the idealized character of our model, we believe that the results are of general 
interest, and some of our predictions could be tested in detail either in extremely clean 
electronic systems subject to strong electric fields or, possibly, in vapors of ultra cold fermonic 
atoms^^. 

This paper is organized as follows. In Section II we introduce the model and set up the 
quantum kinetic equation for the response to an electric field. In Section III we present the 
analytic solution of the kinetic equation and its long-time limit. In Section IV we discuss 
the perturbative regime 7 -C 1. Section V presents the calculation of the transient spin 
Hall current. Section VI describes the inverse (in the sense of Onsager reciprocity) of the 
nonlinear Edelstein effect. Finally, section VH presents a qualitative discussion of the effects 
of disordered and the prospect for experimental observation. 

II. MODEL AND KINETIC EQUATION 

Our model (two-dimensional electron gas with Rashba spin-orbit coupling a) is described 
by the Hamiltonian 

^(P, ^) = 7 ^ [P + eA{t)f - a([p + eA(t)] x cr) ■ z , (1) 

2m 

where 

A{t) = —eEt-k ( 2 ) 

represents an electric field in the x direction, and t is time. The system is assumed to be 
initially (i.e., at time f = 0 ) in equilibrium with a density matrix 

P(P,0) = f+{p)l + f-{p){p X cr) ■ z (3) 
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FIG. 1 . Equilibrium distribution of electrons in momentum space. The inner circle, of radius pF-\- 
has doubly occupied states. The annulus pF+ < pf < Pf- (shaded area) has singly occupied 
states, with spins oriented along the Rashba effective field. The density of electrons in the annulus 
is ~ 2najvF, where n is the total density: only these electrons respond to the application of the 
electric field. 

where 1 is the 2x2 identity matrix, 

9{pf+ -p) ±0{pf- -p) 

2 

and pf+,Pf- {pf+ < Pf-) are the Fermi momenta in the two chirality bands (see Figure 1 ) 
(we work at zero temperature for simplicity) For t > 0 the evolution of the density matrix 
is obtained by solving the equation of motion 

+f[Ff(p,t),p(p,t)] = 0 (5) 

We immediately observe that for p < pf+, i-e. in the region where momentum states 
are “doubly occupied” by electrons of opposite chiralities, one has f-{p) = 0 and hence 
the initial distribution function is proportional to the identity matrix. This commutes with 
the hamiltonian and therefore remains constant in time. This means that the states with 
p < Pf+ constitute an inert background. We will, from now on, focus exclusively on electrons 
in the annulus pf+ < p < Pf- where the states are singly occupied. In this region the initial 
density matrix describes a pure state in which the spin is aligned parallel to the Rashba held 

Bk(p) = apz X p (6) 

The subsequent evolution of this state is caused by the action of a time-dependent “Edelstein 
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field”, which does not depend on p: 


^E(t) = az X eA{t) = —a{eEt)y. ( 7 ) 

The problem is now reduced to calculating the evolution of a single spin, initially aligned 
along B^(p) under the action of the total Zeeman held 

Btoi(p, t) = Br(p) + BE(t) (8) 

This problem is solved analytically in the next section. Once the spin dynamics is solved, 
we can also calculate the charge current response, which is the sum of a linearly varying 
“diamagnetic” term due to the vector potential (this is simply the ballistic acceleration of 
the electrons) and the spin term aS^. 

III. ANALYTIC SOLUTION OF THE MODEL 

Our system of coordinates is shown in Fig. 2 . The hxed direction of the Edelstein held (y) 
is taken as our X3 axis. The original —x and z directions become our xi and X2 respectively. 
We denote by 6*p the angle between p and the standard x axis. Thus the Rashba held forms 
an angle 0 p with the X3 axis and we can write 

Bji(p) = ap(cos0pX3 + sin^pXi), 

BB(t) = —aeEtx 3 , 

Btot(p,i) = (ap cos 0 p — aeEt)x3 + apsin6'pXi, ( 9 ) 

We use the projections of the spin along the X3 axis as the basis for our representation of 
the spin. Then the initial state of the spin of the electron with momentum p (in the range 
Pf+ <p < Pf-) is 

l^p(O)) =cos^| t)+sin^| ;). (10) 

The time-dependent hamiltonian is 

( —ap cos 9 jj + aeEt —an sin \ 

^ ^ , (11) 

—apsinOp +ap cos 6 p — aeEt I 

which is recognized to be the canonical Landau-Zener Hamiltonian^®’^^, which is often used 
to describe the transition probability between two energy levels {—ap cos 9 p -|- aeEt and 
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FIG. 2 . Coordinate system used in the calculations. The spin is initially oriented along the Rashba 
field Br(p). The angle Op between the Rashba field and the X3 axis is the same as the angle 
between p and the x axis. B£;(t) is the Edelstein field, which grows linearly in time, starting from 
0 at the initial time. 

ap cos dp — aeEt in this case), which anticross as a function of time (this problem was hrst 
studied in the context of spin physics by E. Majorana^^. For a pedagogical discussion of the 
model see Ref. 24.) The dynamics depends crucially on whether the rate of variation of the 
energy levels as they cross, eE/{pcos9p), is small or large compared to the magnitude of 
the matrix element that couples the two levels, apsin^p. In the former case the process is 
adiabatic and the spin follows faithfully the magnetic held; in the latter the spin has no time 
to respond to the rapidly changing conditions and remains close to its initial orientation. 
An overall measure of non-adiabaticity is therefore given by the ratio 

_ eE eELg 

ap% Ep 

where Lg = l/{2ma) is the spin precession length (here and in the following we set h = 1). 
The adiabatic regime is characterized by 7 -C 1 and the non-adiabatic one by 7 S> 1. We 
stress that 7 is only an average indicator of adiabaticity: states with = 0 will never be 
adiabatic, no matter how small 7 is, because the matrix element coupling the two levels 
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vanishes when 0p = 0. 

In order to simplify the calculations that follow we express time in units of {apF^/l) 

T 


- 1 . 


t = 


apFy/l 

Then the Hamiltonian (also expressed in units of apF\/^) takes the form 


(13) 


= 


r - Tn 


-A, 


-A. 


r - Tn 


(14) 


where we have introduced the notation 

cos 6-c 

T-p = 


. sm 
A = 


(15) 


\/7 

Notice that Tp is the (dimensionless) time for which the gap between the levels would close 
in the absence of the coupling Ap. The absolute value of Ap is the “residual gap” at the anti 
crossing point. Because the hamiltonian depends on time only via the combination t — Tp 
it is evident that we can represent the solution of the time-dependent Schrodinger equation 
in the form 

iV’pW) = Mp(r - rp)\ t) + Vp{T - Tp)! ;) , (16) 

where the amplitudes Upi^r) and Vpi^r) satisfy the system of equations 

iiipir) = TUpir) - Apnp(r) 

^hp(r) = -rnp(r) - ApMp(r). (17) 

where the dot denotes the derivative with respect to r and the initial conditions are 


Upi-'Tp) = cos(6'p/2), ^'p(-'^p) = sin(6'p/2). (18) 

Two mutually orthogonal solutions of Eqs. (17) are readily found in terms of parabolic 
cylinder functions D{v,z) (see Appendix for details of the derivation) as follows: 

^^(r) = T)(-fA2/2,e-*T\/2r)e-"^p/® 

'^pH'^) = e*^(Ap/A/2)Zl(—iAp/2 — 1, , (19) 


and 


vf{T) = 


( 20 ) 





The constant e has been chosen so that both solntions satisfy the normalization con¬ 

dition 


\up{r)\^ + \vp{T)\‘^ = 1 


( 21 ) 


as they shonld. It tnrns ont that these solntions are precisely the solntions of the classic 
Landan-Zener problem, in which the system is prepared at r = — oo in one the two eigen¬ 
states I t) or I 4 ,). Indeed, making nse of the asymptotic behavior of the parabolic cylinder 
fnnctions^®’^® we hnd 


^ g-i[r2/2+A2 ln(V2r)/2] 


e-^A 2/2 


(r —)■ —cxd) 
(r —)■ -|-oo) 


and 


yW(T) = g*h^/2+A2 ln(^T)/2+7r/4+argr(l-iA2/2)] 


0 (t" —00) 

\/l — e“^^p (r —)■ -|-cx)) 
where r(^) is the gamma fnnction. This leads ns to the identification of 


( 22 ) 


(23) 


pLZ ^ g-^A2 ^ ^-Isin^ep 


(24) 


as the snrvival probability of the initial level, i.e., the probability of staying on the diabatic 
track at the anticrossing point). 1 — is the Landan-Zener transition probability, i.e. the 
probability of avoiding the crossing. From the dehnition (15) of Ap we see that tends 

to 0 for 7 —0 and to 1 for 7 —>■ 00, nnless = 0, in which case it is 1 for all valnes of 7 

(this indicates the complete breakdown of the adiabatic approximation). 

Finally, we seek onr solntion as a linear combination of the two orthogonal solntions (19) 
and (20) that satisfy the initial conditions of onr problem, as stated in Eq. (18). We set 

Up{T - Tp) = ApMW(r - Tp) + - Tp) 

Vp{T - Tp) = ZlpnW(r - Tp) + - Tp), (25) 

and determine A and B from the conditions 

cos ^ = ZlpMW(-rp) + 5pMj,2)(_T-p) 

sin y = ZlpnJ,^n-T-p) + B^v^^\-t^) . (26) 

On acconnt of the relation (20) between the (1) and (2) solntions we easily hnd 

Ap = cos ^[mW(-Tp)]* + sin ^[^^(-rp)]* 

Bp = sin |<)(-rp) - cos ^<)(-rp). (27) 
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FIG. 3. Plot of the Time evolution of the most responsive spin (0p = 7 r/ 2 ) vs r for three different 
values of 7 . Apart from stronger oscillations, which are washed out by angular integration, these 
plots are qualitatively similar to the plots for the time evolution of the total spin polarization in 
Fig. 4. 

This, combined with the expressions of Eq. (19) completes the analytical solution of our 
model. 

In terms of the amplitude Mp(r — Tp) we express the ^/-component of the spin as 

\ (i«p(^ - ^p)r - - ^p)p) 

= |Mp(T-rp)|2 - i. (28) 

At last, the total Edelstein spin polarization is the number of states satisfying pf+ < p < pf- 
(i.e., 2an/vF) times the angular average of 5 'p^ 3 (r): 

■ 1 ) ’ 

where the magnitude of p is approximated as pf- 

Figs. 3 and 4 present plots of the exactly calculated analytical solution as a function of 
time for three representative values of 7: 7 = 0.1 (quasi-adiabatic regime), 7 = 1 (interme¬ 
diate regime), and 7 = 10 (sudden switch-on, or “anti-adiabatic” regime). Fig. 3 shows the 
time evolution of the most significant (i.e., most responsive) spins with 6p = tt/2. Fig. 4 
shows the total spin response integrated over the Fermi surface. 
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FIG. 4. Plot of the Edelstein spin polarization Sy/n in units of a/vp vs r ( 7 r = —v{t)/vF) for 
three different values of 7 . Notice that long-time limits are given by the asymptotic formula (30) 
and are equal to -0.449323, -0.286822, -0.0982924 for 7 = 0 .1,1, and 10 respectively. 


Let us now consider the long-time limit of the solution. Substituting the asymptotic 
forms (22) and (23) into the expression (28) for the Edelstein spin polarization we obtain 


Mp(cX))P 


2 / 771 ^ 7 ^ Se 17B;e‘ 


f+argr 




(30) 


This remarkable result tells us that the Edelstein spin Sy^-p tends, for large time, to a 
constant limiting value, |Mp(oo)p —1/2. This is expected, because in this limit the Edelstein 
held is much larger than the original Rashba held, and the projection of the spin along 
its direction becomes essentially a constant of the motion. However, the limiting value is 
strongly dependent on 7 . For 7^1 the evolution of the spin is generally adiabatic, with 
the exception of states with 6 'p ~ 0 for which sin^^p/y < 1 (see discussion in the next 
section). Thus, with the exception of Op in the immediate vicinity of 0, the spin follows 
the orientation of the total ehective magnetic held, settling in a state with Sy^p = — 1/2 for 
r —)■ 00 . Mathematically, this corresponds to the fact that |Mp(oo)p —)■ 0 for 7 —)■ 0. 

In the opposite limit of 7 S> 1 the evolution of the spin is strongly non-adiabatic. Basically 
the projection of the spin along the direction of the Edelstein held does not have enough 
time to change: it remains equal to the initial value in the limit of 7 —)■ 00 . Mathematically 
this is expressed by the fact that |Mp(oo)p —)■ cos^ ^ for 7 —)■ 00 . Fig. 5 shows the inhnite 
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FIG. 5. Plot of Sy^p{oo) vs. 7 for various values of 6p ranging from vr/S (top dashed line) to 7tt/8 
(bottom dashed line) in steps of tt/S. In all cases, the limiting value at 7 = 0 is Sy^p{oo) = —1/2, 
signifying perfectly adiabatic response. The limiting value for 7 —?• 00 is cos0p, signifying that 
<S'y,p(oo) = 5j;_p(0). The solid line is the angular average of Sy^p{oo). Notice that the transition 
from the adiabatic to the non adiabatic regime occurs at lower values of 7 when Op is close to 0 
and TT (top and bottom lines). In fact, the adiabatic regime is never attained at 0p = 0 or tt. 

time limit of S'j^^p(oo) = |Mp(cx))p — | as a function of 7 . 


IV. LIMIT OF 7 <1 

Let us examine more closely the important limit of 7 -C 1, i.e., weak electric field. For 
a given angle, ^p, the parameter that controls the “adiabaticity” of the dynamics is the 
ratio of the fractional rate of change of the effective field |Bp|/|Bp| to the energy difference 
between the two opposite orientations of the spin in the total Zeeman field, |Bp|. This gives 

I Bp I aeE 

{apF cos 6p — aeEty + {apF sinOpy 

= _2__ ('3X) 

(cos 9p - + sin^ 9p 

The adiabatic regime occurs when rjp -C 1. For small 7 , this will always be the case for 
the states with tt/ 2 < 9p < tt for in this case the denominator is always larger than 1 . 
On the other hand, for states with 0 < ^p < 7 r/ 2 , the denominator reaches the minimum 
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value sin^ when r = cos 6 *p /^/7 = Xp. The condition of adiabaticity is satished only for 
7 -C sin^ 9p or, for a given 7 , when 

I sin 6 'p| 3> . (32) 


Thus, the adiabatic approximation always fails for 6p close to 0, no matter how small 7 is. 
The non-adiabaticity “kicks in” at r = Xp, i.e., at the crossing of the levels: this occurs when 
the velocity of the electrons equals the Fermi velocity. 

These qualitative considerations are conhrmed by an explicit calculation of the adiabatic 
spin response to the Edelstein held. We have 



1 

2 



X 


1/2 


(33) 


and therefore the ^/-component of the spin is given, according to Eq. (28) by 


= 2 - 


y/yr - cos 


(34) 


-^1 + 7 x 2 _ 2y/^T cos 6 'p 

where we have made use of Eqs. (15) for Xp and Ap. The Edelstein spin density is then 
obtained by performing an elementary integration over the angle 6p: 



an 1 1 cos 6 — 

— — , 

vf 2 27r Jo + (a/7t)^ - 2^xcos0 


This integral is elementary and gives 



an 1 

vp 2 


(y7x + 1)E 


477 T 


(vb^+i)^ 


+ ( 7 X - l)K 




(Vh-r+l)^ 


TT 




(35) 


(36) 


where E and K are the standard elliptic integrals^®. 

Fig. 6 shows the Edelstein spin polarization as a function of ^/yx = |u(x)|/uf, where 
v{t) = — —T— is the velocity of the freely accelerating electrons at time t = —. In the 

^ ' m y/yapp j j o yfEjocpp 

linear response regime \v{t)\/vF 1 this formula reduces to 


Sy{t) ~ 


an eEt 
2vf mvF 


2 


aeEt , 


(37) 


where Nq = n/ep = nr/(27r) is the density of states (per spin). This is the standard formula 
for the linear Edelstein effect. Notice that the shortcomings of the adiabatic approximation 
do not show up in this regime, because one is never close to the Landau-Zener anticrossing. 
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FIG. 6 . Plot of the Edelstein spin polarization Syjn in units of ajvp vs r for 7 = 0.1, taken 
as representative of the “adiabatic regime”. The dashed line is the adiabatic approximation of 
Eq. (36). The solid line is the exact solution. The discrepancy between the two curves for yr > 1 
is attributed to the persistent non-adiabatic response of spins with 0p ~ 0 and vr. The dash-dotted 
line is calculated by the improved formula (38) and it is seen to be in much better agreement with 
the exact results. 


When i/qr = 1, (meaning that |n(r)| = vp) a non-analyticity (logarithmically inhnite 
derivative) is present and clearly visible in the plot of Fig. 7. In the next section we show 
that this leads to an unphysical divergence of the spin current. These are all artifacts of 
the adiabatic approximation, and can be cured in an approximate but very effective manner 
by multiplying the integrand of Eq. (35) by the probability 1 — of the Landau-Zener 
transition, i.e., the probability of actually staying on the adiabatic track. The resulting 
formula. 


Sy{.r) 


cm 1 1 cos 6 — 

-/ dO — , : 

Vp 2 27r Jo y/f -k (^r)2 - 2^tcos9 




(38) 


is numerically evaluated and plotted in Fig. 6. This formula is free of pathological behaviors 
at the anticrossing point, since the contribution from the “unresponsive” states with sin6*p < 
^/7 has been suppressed by the Landau-Zener transition probability. 
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FIG. 7. Dashed line: plot of the transient spin Hall current jy oc dSy/dr in units of —eili/( 47 r) vs 
r calculated according to the adiabatic formula (36) for 7 = 0.1: the unphysical divergence occurs 
when ^/yr = 1, that is to say when the drift velocity of the electrons equals the Fermi velocity. 
The solid line is calculated according to the improved formula (38) and it is free of singularities. 


V. SPIN HALL CURRENT 


In this section we calculate the transient spin Hall current, which accompanies the electric 
current. From the Heisenberg equation of motion 

( 39 ) 


we immediately obtain 

= -o^Py^z = -2maJy (40) 

where Jy = is the operator of the spin current. The total spin current density is 

therefore proportional to the time derivative of the total Edelstein spin density. The spin 
current can be calculated analytically, according to the formula 


1 dSy{t) 
2ma dt 



(41) 


In the linear response regime, making use of Eq. (37) we recover the well-known result 


r 



(42) 


Beyond the linear response regime (but still in the quasi-adiabatic regime 7 -C 1) we can 
calculate jy{t) from the time derivative of Eq. (38), or from the exact formulas. The results 
are plotted in Fig. 7. 
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VI. INVERSE EDELSTEIN EFFECT 


The calculations we have done in the previous sections can be straightforwardly adapted 
to the inverse Edelstein effect. This is the reciprocal, in the Onsager sense, of the direct 
effect. We apply a magnetic held that couples to the y component of the spin and varies 
linearly in time: 

By{t)=Tt, (43) 

and we calculate the charge current 


Jx e 



^^n,y 


(44) 


that hows in response. The sum runs over the particles, labelled by n. Since the momentum 
distribution is not ahected by By{t) we immediately conclude that the expectation value of 
vanishes and we are left with 

m 


Jxit) = 2aeSy(t). 


(45) 


Observe that By{t)/2 corresponds to —aeAxit) of our previous calculation, and the “spin 
injection held” r/2 corresponds to aeE. Thus, if we represent the result of our calculation 
for the direct Edelstein ehect in the form 


where we have used the fact that 7 = eE/{ap^p) and = eEt/pp and E is the appropriate 
function of the two arguments, we immediately conclude that the current generated by the 
inverse Edelstein ehect is 


Mt) = 2aeF . (47) 

In particular, in the limit of large times the current tends to a limiting value proportional 
to 

It is interesting to observe that the direct and inverse “conductivities” are 


dSyjt) ^ f ^ ^ ^p{2) (^ 

dE ap% \apy Pf ) Pf Pf J 

dJxit) ^ ( r rt \ et ^(2) / r rt \ 

dr aFp 2app) pp \2a^Pp’2app J ’ 


(48) 


(49) 
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where and denote the partial derivatives of F with respect to the hrst and the 
second argnment respectively. This generalized Onsager reciprocity relation remains valid 
well beyond the linear response regime. The condnctivities of the direct and inverse processes 
are identical if and only if they are evalnated at helds E and T that satisfy the reciprocity 
condition 2aeE = T. This condition is of conrse satisfied in zero field, where onr relation 
rednces to the standard reciprocity relation of linear response theory. 


VII. EFFECT OF DISORDER AND PROSPECTS FOR OBSERVATION 

Up to this point we have completely neglected the effect of impnrity scattering. This 
has allowed ns to obtain an exact and completely analytical solntion. However, it raises 
qnestions abont the possibility of observing the nonlinear effect in realistic system. Cold 
trapped atoms, being intrinsically free of disorder, conld provide an opportnnity to do this. 
The challenge is to find a way to create an artificial Rashba spin-orbit field for cold atoms. 
So far only pnre gange fields (e.g., the eqnal weight combination of Rashba and Dresselhans 
helds) have been snccessfnlly synthesized by exposing the atoms to mnltiple laser helds which 
induces a qnantnm coherence between two hyperhne levels of the atom (the “spin” degree 
of freedom).However, there seems to be no obstacle, in principle, to the realization of 
an artihcial two-dimensional Rashba held, and, in fact, theoretical proposals to this ehect 
have already been pnt forward. 

Let ns fnrther consider the case of electrons in clean systems. Dne to the nnavoidable 
presence of impnrities momentnm is not conserved and the distribntion fnnction is no longer 
constant in momentnm space. In the relaxation-time approximation it shifts along the 
direction of the electric held by a time dependent qnantity mv{t) which eventnally satnrates 
to the Drnde valne mv{t) = —eEr, where r is the electron-impnrity scattering time. A 
reasonable approximation is mv{t) = —eFr(l — which prodnces an Edelstein held 

BE{t) = -aeFr(l - e"*/"). (50) 

The problem is now to calcnlate the spin dynamics of the electrons in this time-dependent 
held, which is no longer linear. If r is snfhciently long the resnlts will be indistingnishable 
(for t <C t) from those obtained in the previons section. The non linearity will be observable 
if the terminal velocity of the electrons is comparable to the Fermi velocity. Alternatively, 
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one could have a very large electric field acting on a system with a not-so-large r. 


The inverse Edelstein effect is more delicate. The charge current will now have contribu¬ 
tions not only from the injected spin {aSy), but also from p^. The latter arises because the 
applied field By changes the distribution of the electrons in momentum space. In fact, under 
equilibrium conditions the contribution would exactly cancel the aSy contribution. One 
way to calculate the effect is to solve the spin dynamics in the presence of a linearly growing 
field By = TtQ(t) in the absence of impurities (which gives us the already calculated current 
aSy) and then take into account the impurities by subtracting the Px current generated by 
the shift in the momentum distribution. For the latter, in the spirit of the relaxation time 
approximation, we assume that it is the equilibrium distribution in a “retarded” magnetic 
field r[t — r(l — 6“*/”^)]. This reduces to the clean result in the limit r —)■ cx) (since the 
external field vanishes for negative times and the equilibrium distribution carries then no px 
current). Whereas in the steady-state regime t 3> r it yields a result proportional to r as 
expected from Onsager reciprocity. Once again, we conclude that the nonlinear effect can 
be observed if the system is sufficiently clean. 


As a final point we wish to comment on what happens in the case that the electrons are 
in a Bloch band with periodic dispersion. The Bloch wave vector k is still a constant of the 
motion. The Edelstein field oscillates in time at the Bloch frequency ojb = where a 
is the lattice constant. This will induce oscillations in both the Edelstein spin polarization 
and the spin current. Depending on whether tag is small or large relative to app we will 
have adiabatic or non-adiabatic response. 
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Appendix: Solution of Eqs. (17) 


To see that the functions dehned by Eqs. (19) are the solution to Eqs. (17) we introduce 
a rescaled time variable z = \/2e“*Tr and rewrite Eqs. (17) as follows 

^ / \ / \ 7 7L r) / \ 

Next, we dehne a parameter z/ = —iAp/2 and transform these equations to the form 

d / \ ^ N Vp(z) 

;r“p ^ + b“p ^ 

dz 2 v— 


(A.l) 


d Vp{z) zvp{z) 


+ Up{z) = 0. 


(A.2) 


dz \/—v 2 ^—v 

By comparing Eqs. (A.2) with the following recursion relations for the parabolic cylinder 
functions D{y,z) [see, for example. Refs. 25 and 26] 


^T>(z/, z) + |t>(z/, z) - vD{v -l,z)= 0, 
^D{u -l,z)- ^D{iy -l,z) + T>(z/, z) = 0, 


(A.3) 


we identify the functions Up{z) = Diy^z) and Vp{z) = ^/—uD{v — l,z) as the solutions 
to Eqs. (A.l) and therefore to Eqs. (17). By returning to the original time variable r we 
recover Eqs. (19) up to the normalization factor. 

To construct the second linear independent solution we notice the following property of 
Eqs. (17). If a pair {mp(t), np(r)} is a solution to Eqs. (17), then {—'Pp(r),Mp(r)} is also 
a solution. Moreover these two solution are orthogonal to each other at any r. Using this 
property one obtains the solution of Eqs. (20) from Eqs. (19). 
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